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LAPLACE EIGENFUNCTIONS AND DAMPED WAVE 
EQUATION II: PRODUCT MANIFOLDS. 

by 

N.Burq & C.Zuily 


Abstract. — The purpose of this article is to study possible concentrations of eigenfunc¬ 
tions of Laplace operators (or more generally quasi-modes) on product manifolds. We show 
that the approach of the first author and Zworski [10, 11] applies (modulo rescalling) and 
deduce new stabilization results for weakly damped wave equations which extend to product 
manifolds previous results by Leautaud-Lerner [ 12 ] obtained for products of tori. 

Resume. — Dans cet article, on etudie les concentrations possibles des fonctions propres du 
Laplacien (ou plus generalement de quasi-modes) sur des varietes produit. On demontre que 
I’approche du premier auteur avec M. Zworski [ 10 , 11 ] s’applique (modulo un changement 
d’echelle) et on en deduit de nouveaux resultats de stabilization pour I’equation des ondes 
faiblement amortie qui generalisent au cas des varietes produits des resultats anterieurs de 
Leautaud-Lerner [ 12 ] obtenus dans le cas de produits de tores. 


1. Notations and main results 

In this work we continue our investigation [9] of concentration properties of eigenfunc¬ 
tions (or more generally quasimodes) of the Laplace-Beltrami operator on submanifolds 
and we study here the very particular setting of product manifolds. 

Let {Mj,gj),j = 1,2 be two compact manifolds. We denote by (M = Mi x M 2 , 5 = 
gi ® 52 ) the product, and by dj (resp. d) the geodesic distance in Mj (resp. M). Let 
go G M 2 and 

S = Ml X {go}. 

For /3 > 0 we introduce 

(1.1) Np = {m= {p, q) e M : dim, S) < /3} = Mi x (g G M 2 : d 2 {q, go) < /3}. 

Our first result is the following. 

Theorem 1.1. — For any <5 > 0, there exists (7 > 0, ho > 0 such that for every 0 < h < 
ho and every solution -0 G of the equation 

{h^Ag + l)i; = F 
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we have the estimate 

(1-2) IIV’llL2(Ar^^) < C{\\i’\\L^{N^^^s\N^^s) + 

As an application of Theorem 1 . 1 , we consider weakly damped wave equations on a 
compact Riemaniann manifold {Ai,g), 

(1.3) {d'f - Ag + b{m)dt)u = 0, {u,dtu)\t=o = iuo,ui) £ x 

where 0 <b £ L°°{A4), for which the energy 

E{u){t) = / {gp{\/gu{t,m),\/gu{t,m))+ \dtu{t, 171)1“^) dvg{m) 

Jm 

is decaying since j^E{u){t) = - jj^h{rn)\dtu{t,m)\^dvg{m) < 0. Let 

uj = U{C/ open : ess inf 6 > 0} 

be the domain where effective damping occurs. We denote by 

QC = {p £ S*M : 3s G M; ^{s)p = G S*u}, 

the (open) set of geometrically controlled points (here d>(s) is the bicharacteristic flow). 
Let 

(1.4) T=S*M\gC, T = n^T 

where T is the trapped set and Ha; the projection on the base manifold A4. 

Our second result is the following. 

Theorem 1.2. — Assume that 

1. there exists a neighborhood V ofT in Ad, a compact Lipschitz Riemannian manifold 
[Mi,gi) of dimension k and a Lipschitz isometry 

Q :V ^ (Ml X R(0, l),g = gi® 52) 

where i3(0,1) is the unit ball in endowed with the (Lipschitz) metric g 2 , 

2. there exists 7 > 0 , c, C > 0 such that 

(1.5) c\z\^'^ <b{Q~^{p,z)) < C\z^'^, V(p, z) G Ml X i?(0,1). 

Then there exists C > 0 such that for any {uo,ui) £ H^(Ad) x the solution u 

to (1.3) satisfies 

E{u)^^‘^{t) < —^(||^io||H2(x) + Ikill^bX))- 

Remark 1.3. — A simpler (but weaker) statement would be to assume 

(i) (Af, 5 ) = (A/i X A/ 2,51 ( 8 ) 52 ), go G A/ 2 , r = S = Mix{5o}, 

(m) cd{m, S)^'^ < b{m) < Cd{m, S)^'^ Vm G Mi x Lf. 

It is classical that for non trivial dampings b > 0, the energy of solution to (1.3) converge 
to 0 as f tend to infinity. The rate of decay is uniform (and hence exponential) in energy 
space if and only if the geometric control condition [2, 7] is satisfied. In [9], we explored 
the question when some trajectories are trapped and exhibited decay rates (assuming more 
regularity on the initial data). This latter question was previously studied in a general 
setting in [13] and on tori in [ 8 , 14, 1] (see also [10, 11]) and more recently by Leautaud- 
Lerner [12]. The geometric assumptions in [9] are much more general than in [12] which 
is essentially restricted to the case of product of flat tori. On the other hand, due to 
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this more favorable geometry, the decay rate in [12] is better than in [9]. Theorem 1.2 
shows that Leautaud-Lerner’s result (the better decay rate) extends straightforwardly to 
the case of product manifolds (Mi x M 2,5 = 51 < 8 ) 52 )- 

Remark 1.4- — 1- According to Theorem 1.6 in [12] the rate of decay in t obtained 

in Theorem 1.2 above is optimal in general. 

2. Theorem 1.1 is a propagation result in the ^;-variable in B{0, 1), and since z is actually 
very close to 0 , the relevant object is 32 ( 0 ) (constant coefficients) rather than 52 ( 2 ). 
The smaller 6, the further we need to propagate in the z variable and hence we better 
the quasi modes we need to consider (due to the worse error factor 

3. The case 5 = 1/2 in Theorem 1.1 is a particular case of our results in [9] (which 

are actually much more general and hold without the ’’product” assumption on the 
geometry). On the other hand, the results in [9] are local, while for 6 < 1/2, the 
estimate (1.2) is non local. Indeed, trying to replace ■0 by x'4’ ^-dd to the r.h.s. 
a term [[h‘^A,x]'ip) which is clearly bounded in by 0{h), giving an error of or¬ 
der S> 1 to the final result. On the other hand, as soon as <5 < 1/2, 

estimate ( 1 . 2 ) is false without the product structure assumption as can be easily 
seen on spheres by considering the eigenfunctions = {xi + 1 x 2 )^ with eigenvalues 
An = n{n + d — 1) = h~‘^ which concentrate in an /ly^-neighborhood of the equator 

E = {x & ; |x| = 1, X 3 = • • • = Xd+l = 0}. 

In this case, we get (/i^A -|- l)en = 0, but 

_ u2S-l 

l|en||l,2(7v^^ ) ~ Ci/ln^^ , ||en||i,2(7V2^^ ) < 026 , n -|-00, 

contradicting ( 1 . 2 ) since 26 — 1 < 0. 

4. No smoothness is assumed on the function b £ Notice however that (con- 

trarily to the results in [9]) the lower bound in (1.5) is not sufficient (at least with 
our approach) and we do need also the upper bound. 

5. As will appear clearly in the proof, we could assume that T is isometric to finitely 
many product manifolds, with possibly different constants 7 , the hnal decay rate 
being given by the largest 7 . 

The paper is organized as follows. We first show how to deduce from Theorem 1.1 a 
resolvent estimate which according to previous works by Borichev-Tomilov imply Theo¬ 
rem 1.2. Then we prove Theorem 1.1 by elementary scaling and propagation arguments. 


2. Prom concentration to stabilization results (Proof of Theorem 1.2) 

According to the works by Borichev-Tomilov [3], stabilization results for the wave equa¬ 
tion are equivalent to resolvent estimates. As a consequence, to prove Theorem 1.2, it is 
enough to prove (see [12, Proposition 1.5]) 

Proposition 2.1. — We keep the geometric assumptions in Theorem 1.2. Consider for 
h > 0 the operator 

(2.1) Lh = -h^Ag-l + ihb, b£L°^{M). 

Then there exist C > 0, ho > 0 such that for all 0 < h < ho 
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for all ip ^ 

Proof. — We start with a simple a priori estimate. Multiplying both sides of the the 
equation 


(2.2) {—h'^Ag — 1 + ihh)p = f. 

by Tp, integrating by parts on Ai and taking real and imaginary parts gives 

(2.3) hf b{m)\p{m)\'^ dvg{m) <\\p\\L2(^j^^\\f\\L2(^M)^ 

J M 

(2.4) h'^ j^gm{^gp{m),Vgp{m) dvg{m) < \\p\\12(^m) + \\p\\L^M)\\f\\L2{M)■ 

l^ow, in the neighborhood y of T we use our isometry 0 and we set 

(2.5) u{p,z) = p{Q~^ip,z)), b{p,z) = b{e~^{p,z)), f{p,z) = f{e~^{p,z)). 
Then from (2.2) we obtain the equation on Mi x 5(0,1) 

{h^Ag + l)u = ihhu — f. 

We can therefore apply Theorem 1.1 and we obtain 

(2-6) \\u\\]^2(^Mix{\z\<h^}) ^ C!\\u\\]^2(^Mix{h^<\z\<2h^) A ^ “/IIl2(Mi x{|2|<2h'5})■ 

On the other hand, from (2.3), and the lower bound in assumption (1.5), we deduce 

(2-7) \\'^\\L^{Mix{h^<\z\}) — '^\\‘f\\L^(M)\\f\\L'2iM) 

while from the upperbound in assumption (1.5), we get 


( 2 . 8 ) 


Lp{Mix{\z\<2hS}) - 

Mix{\z\<2h^} 


Gathering (2.6), (2.7) and (2.8) we obtain, 

+ Ch^^-^‘ 


(2.9) 


L^{M) 


Optimizing with respect to <5 leads to the choice 26 = which gives 

ll«lli2(MixB{0,l)) ^ Ch~^~^\\p\\L 2 ^M)\\f\\LHM)+Ch~‘^~^\\f\\l 2 (^j^y 
According to (2.5) this implies 

(2.10) ||(/?||i2(y) < C/i ^ ||7 ’IIi,2(_m)II/IIi,2(_a 4) + G/i ^ Il/lli2(x)- 

We can now conclude the proof of Proposition 2.1 by contradiction. If (2.1) were not 
true, then there would exists sequences pn S 5^(A4), fn G L^(A4), 0 < /i„ —>■ 0 such that 

{ — hnAg — 1 + ihnb)pn = fni ll7’n||L2(^) > ^ , y ll/n||L2('JVfP 

h 

Dividing pn by its L^-norm, we deduce 

IH—^ 

(2-11) \\pn\\L‘^[M) — \\fn\\L‘^{M) — ^^^), n — ^ +00, 
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and from (2.10) we get 
(2.12) 


lim ||<y7n||i,2(y) = 0. 

On the other hand, the sequence {}Pn) is bounded in and extracting a subsequence, 

we can assume that it has a semi-classical measure /r (see e.g. [5, Theoreme 2]). We recall 
that it means that for any symbol a € 

^hrn^(a(x,/i„Il3,)(/?„,(/?„)^2(_yv,) = 


Here, since we work locally, we quantize the symbols a G C^{T*W^) by taking first 
(f) G C'^(M'^) equal to 1 near the x-projection of the support of a and 


a{x, hDx)u 


1 

(27r/i)” 



y'^'^a{x,i)(t){y)u{y)dyd^. 


It is classical that modulo 0{h°°) smoothing operators, the operator a{x,hDx) does not 
depend on the choice of (f). 

From (2.4), the sequence (i^n) is oscillating and hence any such semi-classical defect 
measure has total mass 1 = lim„_}.+oo ||¥^n||L 2 (_A 4 ) (see [5, Proposition 4]). From (2.3) and 
(2.11) we also have (notice that \b\ < C\h\^/'^) 

( h A — ihnb^Pn “t“ fn — o(hn)i2^ 


and consequently (see [6, Proposition 4.4]) the measure // is invariant by the bicharacter¬ 
istic flow. Since from (2.3) it is 0 on S*uj, we deduce by propagation that it is also 0 on 
QC^ and hence from (2.12) it is identically null, since S*{M.) = T U QC. This gives the 
contradiction. □ 


3. Concentration properties (Proof of Theorem 1.1) 

Recall that we have {M,g) = (Mi x M 2,51 ® g 2 )- The proof of Theorem 1.1 follows, 
after taking scalar products with Laplace eigenfunctions in Mi, from a rescaling argument 
and standard (non trapping) resolvent estimates in M 2 . When the metric 52 is hat, the 
scaling argument is straightforward, while it requires a little care in the general case (see 
Lemma 3.5). 

Let B{qo,r) C M 2 be the ball (for the metric ^ 2 ) of radius r > 0 centered at qq. 

Proposition 3.1. — For any <5 > 0, there exists C > 0,/iq >0 such that for every 
0 < h < ho, every r G M, every solution U G H‘^{M 2 ), G G L^(M 2 ) of the equation on M 2 

(-A,,-r)C/ = G 

we have the estimate 

(3T) \W\\L'^[B{qQ,h^)) ^ \\G\\l'^ [B{qo,2h^)'^ ■ 

3.1. Proof of Theorem 1.1 assuming Proposition 3.1. — Let (e^) be a sequence 
of eigenfunctions of the Laplace operator on Mi with eigenvalues — forming an L^(Mi) 
orthonormal basis. For ^ G L^(M), we set ipn{Q) = (V’(') 9 )) Then we have 

V'(P) q) = '^niQ)enip) and it is easy to see that with the notations in ( 1 . 1 ), for r > 0 

(^•^) 11''A 11L2 (TV,.) = \\'^\\L'2{MixB{qo,r)) ~ ^ \\'‘l^n\\L2[B{qo,r)- 

neN 
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Now taking the scalar product of the equation (1.2) with e„ we see easily that {—h'^Ag^ + 
h'^Xn ~ l)'0n = Fn which can be rewritten as 

{-Ag^ - T)^n = h~‘^Fn, T = h~‘^ - A^. 

Applying Proposition 3.1 to this equation yields 

Taking the sum in n and using (3.2) we obtain the estimate (1.2). 

3.2. Proof of Proposition 3.1. — Since the problem is local near qo, after diffeomor- 
phism we can work in a neighborhood of the origin in and we may assume that the 
new metric g satisfies g\z=o = Id. Then we make the change of variables z ^ x = ^ and 
we set u{x) = U{h^x), F{x) = G{h^x). We obtain the equation on u 

{-Agh - h^^T)u = h^^F, 

where g^ is the metric obtained by dilatation g^{x) = g{h^x). The family (g^) converges 
in C°° topology to the flat metric go = Id. Proposition 3.1 will follow easily from 

Proposition 3.2. — Consider a family (gn) of metrics on 5(0,2) C which converges 
in Lipschitz topology to the flat metric when n —>■ +oo. Then there exists C > 0, Nq > 0 
such that for every n > No, r G M, u G 5^(5(0,2)),/ G L^(5(0, 2)) solutions of the 
equation on 5(0,2) 

-t)u = f 

we have the estimate 

(3.3) II'*^IIl2(b(o,i)) < C'(||u||j;^2(5(o,2)\b(o,i)) + Yq7pjY7^ll/llL2(s(o,2))) 

(notice that since gn converges to the flat metric the choice of the metric to define the 
52 -norms above is of no importance). 

Remark 3.3. — Proposition 3.2 is standard for the fixed metric go = Id (see e.g. [ 6 , 
Section 3]), as the annulus {x : 1 < |x| < 2} controls geometrically the ball 5(0,1). As a 
consequence, in the special case of [12] when g = go (and hence gn is also the standard flat 
metric), the proof of Theorem 1.1 is completed. In the general case, we only have to verify 
that the usual proof can handle the varying metric through a perturbation argument, 
which is precisely what we do below. It is worth noticing that the proof belows implies 
that the propagation estimates involved in exact controlability results which are known 
to hold for metrics, see [4], are actually stable by small Lipschitz perturbations of the 
metric. 

For r > 0 we shall set Br = 5(0, r) C 

To prove Proposition 3.2 we argue by contradiction. Otherwise, there would exist 
sequences, an —>■ + 00 , (t„) C M, (u„) C 5 ^( 52 ), (/n) C 5 ^( 52 ) such that 

(3.4) ( — Ag^^ — Tn)Un = fn, 

(3.5) 1 = \\Un\\L‘2{Bi) > ^(ll'*^n||l, 2 ( 52 \Bi) + |l/2 IIII(^ 2 )) 

We now distinguish three cases 

— lim infTn = —oo (elliptic case) 

~ {Bn)NeN bounded (low frequency case) 
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— lim sup„^_,_Qo = +oo (hyperbolic case) 

In the first case, working with a subsequence we may assume that lim„,^+oo Tn = —oo. 
Let C € C^{B 2 ) equal to 1 on - 63 / 2 - Multiplying (3.4) by Qun, integrating by parts and 
taking the real part gives 


(3.6) 



iCUn))) 


Crn\Un\‘^) dVg^^ 


< \\Un\\L^{B2)\\fn\\L^(B2) 


which implies (after another integration by parts) 


(3.7) 

Since Ag^ 
Tn -00 


J C5o-„ Vg^^ttn)) (t„C + 2 

< \\Un\\L'^{B2)\\fn\\L'^{B2) = n + 00 . 

is supported in {1 < |x| < 2} and ||ttn||L 2 (i<|a;|< 2 ) = o(l), we deduce if 


lim 

n—>-+oo 


C\Un\‘^dx = 0 , 


which contradicts (3.5). 

In the second case (low frequency), we can assume (after extracting a subsequence) that 
Tn ^ T and (3.7) shows that the sequence (rinlsg/a) is bounded in Hence, (after 

taking a subsequence), we can assume that it converges weakly in H^{B^/ 2 ) (and hence 
strongly in L‘^{B^/ 2 )- Due to the convergence of the family of metrics, we get 


k 

~^gcrn^n — ~ A^Un + o(l)_f/-l j (^0 = ^ ^ dj ), 

i=l 


and according to (3.5) this implies that the limit u satisfies 

(-Ao -t)u = 0 in V'{B^/ 2 ), u |i<|^|<3/2= 0 . 

Uniqueness for solutions of second order elliptic operators implies that rt = 0 which is 
contradictory with the strong convergence of {un) in L'^{B^/ 2 ) and (3.5). 

Finally it remains to study the last case (hyperbolic). Taking a subsequence, we can 
assume —)• +oo. Moreover dividing both members of (3.4) by Tn we see that Un is 

solution of an equation of type {P{x,Tn ^ D^) — l)Mn = Tn^fn ^ 0 in L^(H(0,2)). The 
sequence {un\\x\< 3 / 2 ) has a semi-classical measure u with scale 


(see the end of Section 2 for a few fact about these measures). Notice that this new semi- 
classical parameter hn has no relationship with the parameter h in Theorem 1.1 First 
of all multiplying both sides of (3.7) by hn = and using the fact that ||rtn||i, 2 ( 52 ) is 
uniformly bounded we deduce that there exists C > 0 such that 

(3.8) hn\\'^xUn\\L^{B 3 / 2 ) - Vn € N. 

Using again (3.7) shows that the sequence Un\\x\< 3/2 is /in-oscillatory (and hence the 
measure v has total mass 1 = limn_>.+oo ll^^*lli 2 (B 3 ^ 2 ))- Now setting Dn = det((5(o-„)ij) we 
can write 


(3.9) 








2A 


(a.D„)8y. 

*J =1 
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The only point of importance below will be that 

(3.10) lim llfi'ai -ll^n - 1 ||h^i.°°(B2 ) = 0. 

Proposition 3.4- — The measure v is supported in the set {(x,C) : |C| = 1} 
invariant by the bicharacteristic flow associated to the metric qq: 

2C ■ Va;l^ = 0. 

The contradiction now follows since by (3.5) we have 0 which im¬ 

plies that i^|i.i<|a;|<i .9 = 0 and by propagation that v\\x \<^/2 = 0- It remains to prove 
Proposition 3.4. 

Proof. — We have for a with compact support (in the x variable) in B{0, 2), 

(a(x, hnDx){h.iri^gi^^ ~ ^)Un, Un) j ^2 ~ (1) + (2 -|- (3), 

( 1 ) — (ai^X, hnDxflhj^^o \')Un,Un^ 

(3.11) (^) — ^ y ((fftfi, dij) hndjUn, hndiU (^X, hxiDx'jUnj ^2 ^ 

hj 

(3) — hfi 'y ^ (nri {^9<j„d)iDnhndjUn.i a {x , hnDx'jUn^ j^2 ■ 
v 

On one hand, using the symbolic calculus, the term (1) tends to 

(z.,(|Cp-l)a(x,C)). 

Now using (3.8) and (3.10) we see easily that the terms (2) and (3) tend to zero when 
n —>■ -|-oo. On the other hand, the l.h.side in (3.11) is equal to 

(fli.X, hfiDx) fm '^ri) jj2 • 

and according to (3.5) tends to 0. We deduce 

Vo e {n, (ICP - l)a(x, C)) ^ supp (n) C {(x, C); |C|^ = !}• 

To prove the second part in Proposition 3.4, we shall use the following lemma 

Lemma 3.5. — Let a £ and b £ Then 

(3.12) \\[a{x,hnDx),b]\\c(L2) < Chn\\Vxb\\L^. 

Proof. — The kernel of the operator [a(a;, hnDx), 6] is equal to (here f £ C^(R^) is equal 
to 1 on the x-projection of the support of a) 

K{x,x') = - —L— f '^\{xX){bix)-Kx'))(t>{x')dC„ 

{2-Khnfl 4cgR'= 

which is for \x — x'\ < hn (since the support of a is compact) bounded by 

(3.13) Ch->‘\\S/xb\\L^\x-x'\, 

while for |x — x'| > hn we can integrate by parts using the identity 





1 


which gives 
K{x, x') = 

{2TThnY 

and hence gives the bound for any iV € N, 


{x - x’) • 


e 


n .1 ,19 

i\x — x'r 


N 


a(x, (){b{x) — b{x')4>{x')d(^, 


(3.14) 


^ i^N—k 


\x — X' 


It follows from (3.13) and (3.14) that 

f \K{x,x')\dx+ f \K{x,x')\dx'< Chn\\'^xb\\L°°- 

JR'* JRfc 

Then Lemma 3.5 follows from Schur’s lemma. 

Denoting by [A, B] the commutator of the operators A and B let us set 

C = ( [0(2;, hnDx), {h^A.g^^ — 1 )] Uri: Un) j^2 ■ 

hfi 

Then we can write using (3.4) and (3.9), 

C = J-( a{x,hnDx),hlfn],Un )^2 = (1) + (2) + (3) 
hn 

( 1 ) — ~ ( [®(3J) hfiDx), (/l^ Ag 1 )] jj2 j 

K 

k 

(2) = ^ ^ ( [® (^) Dx), (9j ( (o’! ^ djl)h^diy''^Un)Uri^ jj2i 

3,1=1 

i ^ ~ ~ 1 

(3) = ^ ^ hn(^\ci{x, h^Dx), g^^{djDn)hndi^Un,Un) J^ 2 - 


(3.15) 


h. 


By symbolic calculus, the term (1) is modulo an 0{hn) term equal to 

(Op({a(x,C), \Cf})Un,Un)^2, 

where {,} denotes the Poisson bracket, and hence tends to 

{u, {a{x, C), ICP}) = (2C • a). 

Let us look to (2). Each term in the sum can be bounded by 

1 


h, 


-1 ( |^fl(x, hnDx), hndj^ (.dan ^jl'jhndi'Un , Un) 


L 2 l 


+ ~ (^^Cli^X, hnDx), dji^)hndiUn, hndjUn)\ 


□ 


T~ li^idan : ^nl^x'), hndl)un, hndjUn)\- 

K 

By the semiclassical symbolic calculus and Lemma 3.5 the norms in C{L?‘) of the opera¬ 
tors \a[x,hnDx).,hndj\ and \a{x,hnDx),g^}n ~ ^ji)] bounded respectively by Chn and 
ChnWVxgiiW 1^00 where C is independent of n. Therefore using (3.10) and (3.8) we deduce 
that (2) tends to zero when n goes to -|-oo. 
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Unfolding the commutator and using (3.5), (3.8) we see that the third term in (3.15) is 
a finite sum of terms which are bounded by C\\djDn\\L°° ■ We deduce from (3.10) that (3) 
tends to zero when n goes to +oo. 

Now, opening the commutator we see that the r.h.s. in the first equation in (3.15) is 
equal to 

, hnDx')hjifni^2 {x, hnDx^'^n) ^2- 

hn hn 

These terms are bounded by Chn\\fn\\L' 2 {B 2 )\\'’^n\\L^{B 2 ) tend to zero when goes to 

1 

+00 since, according to (3.5), ||un||L2(B2) uniformly bounded and = o{Tn) = 

o{h~^) . This ends the proof of Lemma 3.4, and hence of Proposition 3.2 □ 
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